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This article reviews classical and quantum aspects of critical phenomena in gravitational 
collapse. We pay special attention to the origin of the scaling law for black hole mass, 
and to phase transitions in which black hole formation turns on at finite mass. We 
present some new results for perfect fluids with pressure proportional to density. 

1 Introduction 

The discovery, by ChoptuikI of critical point behavior in gravitational collapse is 
one of the most significant achievements of numerical relativity to date. Efforts to 
understand the Choptuik phenomena have opened an exciting area of research in 
general relativity, one which lias-benefited from the interplay between mathematical 



The problem addressed by Choptuik was suggested by the work of Christodoulou 
on the gravitational collapse of spherically symmetric configurations of massless, 
minimally coupled scalar fields. Christodoulou proved the global existence and 
uniqueness of regular solutions for initial scalar field configurations in a neighbor- 
hood of trivial initial data (i.e., flat, empty space). When the initial data is strong, 
in a precise sense, he also established that the black hole mass has a positive lower 
boundfl Choptuik examined the intermediate regime numerically. By construct- 
ing interpolating families of solutions depending on a single parameter r\, Choptuik 
demonstrated that there generally exists a critical value 77* such that S[r/ < 77*] are 
solutions in which the scalar field disperses to infinity, and S[tj > 77*] are solutions 
in which the field collapses to form a black hole; he called these subcritical and 
supercritical solutions, respectively. He also observed two important properties of 
near critical evolutions S[?7 ~ ry*]: 

1. For each family of initial data, the masses of black holes formed in marginally 
supercritical collapse obey a scaling relation 
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M BH = K\t)-ri*\ 



(1) 
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where [3 w 0.37 is universal, i.e., independent of the initial data, although the 
constant if is a family dependent. 

2. The evolution of the scalar field is characterized by scale-invariant echoing 
with a period, in logarithmic time, given by A « 3.4. 

These two observations led Choptuik to speculate that precisely critical evolutions 
approach a unique solution in which scale-invariant echoes accumulate at a massless, 
central singularity. 

It is now known that critical point behavior is a general feature of gravitational 
collapse. Consider the evolution of a physical system consisting of matter coupled 
to gravity, or simply of the gravitational field itself. One can think of the equations 
which govern this system as a map £ from the space of initial data into the space 
of solutions. Suppose the strength of the initial data I[rj\ is characterized by the 
single tunable parameter r\ (taken to be in the range < i] < oo for concreteness), 
then one constructs a family of solutions S[rf\ by 

£ : I[rj\ - S[rj\ . (2) 

Provided weak cosmic censorship b holds, critical point behavior may be expected 
whenever these solutions interpolate between black hole formation and stable, reg- 
ular solutions. c Some of the questions which then arise are: 

1. Is cosmic censorship upheld in all families of interpolating solutions? 

2. What is the origin of the scaling relation for the black hole mass? 

3. Does black hole formation always turn on at infinitesimal mass? 

4. Is the critical point always associated with scale invariant solutions? 

5. What are the semi-classical corrections to the critical point behavior? 

In this contribution, we provide only a flavor of the work that has been done to 
answer these, and other questions. Our presentation is not intended to be complete, 
rather it focuses on issues that were discussed in the parallel session at the 8th 
Marcel Grossmann meeting (MG8). A more complete review, including technical 
details of the many interesting results which have been obtained, has recently been 
prepared by Gundlach.@ 

Nevertheless, we would be remiss if we did not mention some of the notable 
contributions which were not represented at MG8. All of the work to date has been 
restricted to spherical symmetry except for the study of axisymmetric collapse of 
pure gravitational waves by Abrahams and EvansB They demonstrated the exis- 
tence of critical point behavior and derived a scaling law for the black-hole mass 
similar to that in Eq. (|]) with (3 w 0.37. Despite the extreme complexity of the 
analysis, they also presented tentative evidence for scale invariant echoing of the 

fey weak cosmic censorship we mean that, in the evolution of generic initial conditions, singu- 
larities are always hidden behind an event horizon. 

Minkowski spacetime is a special case where the material disperses completely 
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metric on the symmetry axis of the spacetime. The black-hole mass scaling expo- 
nent is intriguingly similar to that observed in scalar field collapse, although the 
echoing period, A w 0.6, of the solution is quite different. 

The initial explanation of the scaling relation (|l|) for black hole mass was made 
possible by the work of Evans and Coleman.EJ They showed that the critical solu- 
tion, in gravitational collapse of radiation fluid, is continuously self-similar by ex- 
amining one parameter families of interpolating solutions, and direct construction 
of the intermediate self-similar attractor. They indicated that the critical exponent 
might be derived by considering linearized perturbations about the critical solution. 
This suggestion was followed up by Koike et al. a who argued that the critical expo- 
nent, determined to be (3 ~ 0.36 by Evans and Coleman, is directly related to the 
largest Lyapunov exponent of perturbations about the critical solution. (This idea 
was also explored by Eardley and Hirschmann liil in the context of complex scalar 
field collapse.) In this way Koike et al. provided the first direct computation of the 
critical exponent as (3 — 0.355 for radiation fluid. Maison E_Hl then extended these 
results to more general equations of state for the perfect fluid. 

It is worth summarizing the argument which relates, the Lyapunov exponent 
to the critical exponent; the discussion follows Maison.II 2 ] The general spherically 
symmetric line element can be written as 

ds 2 = -a 2 dt 2 + f~ x dr 2 + r 2 (d9 2 + sin 2 9 d(f> 2 ) . (3) 

The self-similar critical solution has a = a*(£), and / = /*(£), where £ = r/t and 
the time coordinate is normalized so that t = at the singularity in the critical 
solution. An important feature of this solution is the existence of a sonic point £ sp 
where the fluid moves at the speed of sound relative to £ — £ sp . Now, consider the 
linear stability of these solutions against perturbations with the general form 

-(v-rh)M)t- x , (4) 

where A is a positive constant. These perturbations arise as small deviations from 
precisely critical initial data — consequently, they are proportional to (rj — 77,). 
There is a unique value of A such that the perturbations are regular both at r = 
and at £ sp . An apparent horizon in the perturbed spacetime is given by 
such that 

f*(Zh)-(v-v*)A(thK x = o. (5) 

The mass Mbh inside the apparent horizon is related to its radius r^ by Mbh = 
r7j/2. Substituting th = rh/£h into Eq. (||), we can solve for Mbh as 

M B H = r ft /2oc|r ? -?7*| 1 / A . (6) 

The scaling exponent is then read off as (3 = I/A. 

In addition to these works, Gundlach has made significant contributions to the 
analytic understanding of critical phenomena when the critical solutions have self- 
similar echoes. He has termed this symmetry discrete self- similarity (DSS), and 
has constructed solutions with DSS by writing the metric as a Fourier series,|-and 
recasting Einstein's equations as an eigenvalue problem for the echoing period.0 As 
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in the case of perfect fluids, where Maison was able to predict the critical exponent 
for situations which had not been numericalhz.. explored, the power of Gundlach's 
analysis was demonstrated by his prediction £3 that there should be periodic fine 
structure on the mass scaling relations such that 

M BH oc \n-rt, fe*V n ^ I) , (7) 

where 'J has period A/(2/3) ~ 4.6|_in In \v ~ V*\- Similar arguments were presented 
independently by Hod and Piran liJ who confirmed the existence of this effect nu- 
merically. 

In collaboration with Martin-Garcia, Gundlach considered charged scalar field 
collapse and again correctly predicted the scaling relation that has been observed 
for the final black hole charge .□ Recently he has extended his analysis to include 
small deviations from spherical symmetry, and has derived a scaling relation fox the 
angular momentum parameter of spinning black holes near the critical point EjEj 



3 



The remainder of this paper is organized as follows. In subsection 2.1 
present a simple argument for the periodic fine structure discovered by Gundlach 
and independently by Hod and Piran.H Interesting new phenomenology is observed 
when scale invariance of the underlying mathematical equations is broken. This was 
first demonstrated in the work of Choptuik et al. by considering the gravitational 
collapse of a YangpMills field where the SU(2) charge introduces a fundamental scale 
into the problem. Il3 In subsection |2.2| we discuss similar results for massive scalar 
field collapse; in particular, we focus on a simple criterion to determine when the 
mass term is important in critical evolutions. We then outline the results which 
have been obtained by Choptuik et al. whci-have shown that the approximate 
black hole solution discovered by Van Putten EJ is unstable. This result also has 
implications for massless scalar field collapse since van Putten's solution may lie 
at the threshold of black hole formation. Some preliminary results from a study 



of fluid collapse with p — kp are presented in subsection 2.4; in particular, the 
critical exponents are tabulated for several values of k between zero and unity, 
and graphical evidence for self-similarity in near critical collapse with k = 0.95 
is presented in Fig. |^. Concluding the section on classical results, we discuss the 



implications of critical point behavior for cosmic censorship in subsection 2.5. To 
fully understand the physical significance of critical point behavior in gravitational 
collapse requires the inclusion of quantum effects into the picture. Unfortunately, 
we do not have a complete theory of quantum gravity, so the best we can do is 
address some model problems. Several interesting results have been obtained in 
this context. Peleg et alJtB have investigated a 2-dimensional dilaton theory of 
gravity in which they demonstrate the existence of critical behavior and scaling at 
the classical level. When one loop quantum effects are included into the model, 
a mass gap is observed at the threshold of black_hole formation. These results 
are reinforced in similar work by Ayal and Piran,E3 who have considered a semi- 
classical model for spherically symmetric scalar field collapse which is outlined in 



subsection 3.1. General arguments can be brought to bear on the four-dimensional 
problem when the critical solution is continuously self-similar. We present a brief 
summary of a recent analysis by Brady and Ottewill c3 which suggests the existence 
of a quantum mass-gap at the threshold of black hole formation for sufficiently stiff 
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perfect fluids. Finally, we conclude with a brief discussion of future directions. 

2 Classical results 

2. 1 Massless scalar fields 

Since the initial discovery of critical phenomena, spherically symmetric scalar field 
collapse has been one of the most studied systems. In contrast to perfect fluid 
collapse the critical solution is not self-similar in the familiar sense, rather it has an 
infinite train of self-similar echoes which accumulate at the singularity. In precise 
terms, scale invariant quantities, such as / in Eq. (||), satisfy 

/(C > ln|t|) = /(C,ln|t|-nA/2) (8) 

where n is an integer and A « 3.4 is the period computed by Choptuik for the 
scalar field. In words, scale invariant quantities are periodic in In \t\ on surfaces 
of constant £. The discrete self-similarity of the critical solution has implications 
for the mass scaling law in Eq. (Q) — there is periodic -fine structure superimposed 
on the power-law. This was first noticed by GundlachJ13 and subsequently verified 
numerically by Hod and Piran.B 

A simple argument demonstrates this result. The mass scaling law arises be- 
cause a single unstable mode governs the way that marginally supercritical solutions 
run away from criticality. This mode can be written as 

-(v-v*)t- x fi(tMt\) (9) 

where /i(£,ln|t|) inherits the symmetry of the background, namely that it is peri- 
odic in In Let th satisfy 

f*-(r)-V*)th X h = (10) 



for some fixed £ = £h- Substitute th = rhjSjx into Eq (10) and solve for rh to get 



lnr h = - In |r?- 77*| In r-ft) , (11) 

where ty(£h, hir^) is periodic in lnr^ with period A/2. If ^ = we have a precise 
power-law relation, however ^ is generally non-zero so an approximate solution for 
r h is 

\nr h = -ln|ry-?7*| + *(e/ l ,A- 1 ln|rj-7 ? *|) . (12) 



Since ^(^ft, hir^) is periodic in lnr/j, Eq. (12) implies that rh, and equivalcntly 
the black hole mass, has periodic fine structure with period AA/2 in \n\rj — r}*\ 
superimposed on the familiar power-law. 

A natural extension of Choptuik's results which was considered by Gundlach 
and Martin-Garciajj and also by Hod and Piranjiil is scalar electro-dynamics cou- 
pled to gravity. The different scaling of the mass and charge in the self-similar 
echoing regime suggests that the significance and influence of the charge decreases 
during a near-critical evolution. As r\ — 77* — > the black-hole charge tends to zero 
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faster than its mass, i.e., Qbh/Mbh —>■ as r\ — ► rj* where Qbh is the black hole 
charge. Gundlach and Martin-Garcia computed the critical solution and predicted 
that the charge should scale as 

n I*, *, |0.883±0.007 /ifl 

WBH OC \T] — 77*| . (16) 

This scaling relation was confirmed by Hod and Piran in their numerical simulations. 



2.2 Type I phase transitions 

By analogy with statistical mechanics, phase transitions in which black-hole forma- 
tion turns on at infinitesimal mass have been termed Type //transitions. Choptuik, 
Chmaj and Bizon E3 observed that black hole formation sometimes turn on at a 
finite mass in the collapse of Yang-Mills field coupled to gravity. They called this a 
Type I transition since the order parameter, the black hole mass, is not continuous 
at the critical point. The critical solution in this sector of the theory is the Bartnik- 
McKinnon soliton, an unstable static solution of the Einstein- Yang-Mills (EYM) 
equations. The mass gap at the threshold of black hole formation is approximately 
equal to the mass of the Bartnik-McKinnon solution. A fundamental difference 
between the EYM system and the massless scalar field is the presence of a length 
scale, the Yang-Mills charge. 

It has been argued that Typc-II,tuansitions in scalar field collapse should be sta- 
ble to the introduction of a mass.Q'E3E3 Gundlach has also presented a mathematical 
argument to this effect in his recent, review 3 Nevertheless, it was demonstrated by 
Brady, Chambers and Gongalves E-3 that Type I transitions can occur in massive 
scalar field collapse. In contrast to EYM collapse the mass gap at the threshold 
of black hole formation is not universal; marginally supercritical black holes have 
masses ranging from ~ 0.3/i _1 to ~ O.QfX , where fj, is the scalar field mass. A 
plausible explanation for this is provided by studying the critical .solutions in Type 

I transitions; they are oscillating soliton stars. Seidel and SuenEil have discussed 
these solutions in some detail. In particular, they demonstrated that there is a fam- 
ily of solutions parameterized by the mass and the effective radius of the star. A 
schematic representation of the mass-radius curve is shown in Fig. |l]; note that the 
maximum mass of such a star is ~ 0.6fi , and solutions to the left of the maximum 
are unstable. The numerical evidence suggests that all of the unstable solutions can 
act as critical point solutions, perhaps explaining the range of values for the mass 
gap. 

Why do some initial configurations lead to Type I transitions while others lead 
to Type II? A complete answer to this question is unavailable at this time, however 
there is a simple criterion which provides some guidance. Let A denote the radial 
extent of the initial shell of scalar field, i.e., its thickness. Type I transitions occur 
when the radial extent A of the initial pulse is bigger than the Compton wavelength 
of the scalar field, that is 

A/x » 1 . (14) 

This observation, and the evidence in support of it that is presented in Chambers' 
articleJ23 validates local arguments about the relevance of the mass term in Type 

II transitions. 




»~ 

Effective Radius 



Figure 1: A schematic representation of the mass-radius curve for oscillating soliton stars. (See 
Seidel and Suen Cj for a quantitatively correct version.) The maximum mass that one of these 
configurations can have is ~ 0.6fi~ 1 . Solutions on the unstable branch of the mass radius curve 

may act as critical solutions 



It is likely that unstable, confined solutions play a role in critical point behavior 
of other matter models. Indeed, this is expected in the astrophysical context where 
stars which exceed the Chandrasekhar limit must either shed some material or 
collapse to form a black hole. The possible end-points of gravitational collapse in 
this context are: (i) black hole formation, (ii) complete disruption of the star, or 
(iii) a dead star with a mass less than ~ 1.4M (where M Q is the mass of the sun). 
In the most general context, phase transitions between any pairwise combination of 
these types can be expectedj-although all of the possibilities need not occur. Recent 
work by Chmaj and Bizon E3 lends support to this viewpoint. They have shown 
that stable skyrmion solutions are a possible endpoint of gravitational collapse in 
Type I transitions, i.e., phase transitions are observed between stable stars and 
black holes. 

2.3 Brans-Dicke Theory 

Critical point behavior in Brans-Dicke theory has been studied by Liebling and 
Choptuik.EEl The critical solution exhibits continuous self-similarity or discrete self- 
similarity depending on the strength of the coupling of the Brans-Dicke field to the 
matter (taken to be massless scalar field by Liebling and Choptuik). In a broader 
context, van Putten has presented a one-parameter family of static solutions to 
vacuum Brans-Dicke theory which he suggests might be thought of as approximate 
black holes L3 These solutions can be arbitrarily close to the exterior Schwarzschild 
spacetime, but are globally regular. An important application of solutions like this 
could be to provide approximate boundary conditions at the surface of a black hole. 
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Choptuik, Hirschmann and Liebling have argued that, in addition to the static 
solution behaving like a black hole, the solution should also be stable. Unfortunately, 
a linear stability analysis indicates that van Putten's solutions have an unstable 
mode. Using the code developed previously, Choptuik et al. also studied the non- 
linear stability of these solutions. They evolved initial data corresponding to the 
static solution and small additive perturbations. Generically, the solution either 
collapses to a black hole, or disperses depending on the sign of the perturbation. An 
interesting point was raised by van Putten during the session on critical phenomena 
at MG8. He noticed, in animations of the perturbed solutions, that the instability 
is evident only after the perturbation reaches the center of symmetry and reflects 
back out toward larger radii. One may therefore wonder if the solutions might be 
stable when absorbing boundary conditions would be applied at the origin; if so, 
the solutions would become interesting as approximate black holes once again. 

Choptuik et al. have shown that van Putten's solutions have a single unstable 
mode, and they either disperse or collapse when perturbed-.-_.This suggests that 
the solutions may be black-hole threshold critical solutions It is intriguing to 
consider this possibility, since it would imply that there is a basin of attraction in 
the space of initial data, for massless scalar field collapse, in which van Putten's 
solutions represent the critical solution. Once again this begs the question: what 
makes one solution a critical attractor, and another not? Clearly, there is a gap in 
our understanding of critical phenomena. 

2-4 Perfect fluid collapse 

Despite the early study of critical phenomena in the collapse of radiation fluid by 
Evans and Coleman,-- no results have been available on the evolution of perfect 
fluids with the general equation of state p — kp, where < k < 1 is a constant, p is 
the energy density of the fluid, and p is the pressure. Solutions expected to be at 
the threshold of black hole formation have been computed by Maison __ and Koike 
et alB when < k <J 0.899. In each case a single unstable mode has been found, 
and a scaling exponent for the black hole mass in slightly super-critical collapse has 
been predicted. 

Brady and Cai __ have developed spherically symmetric code to study this 
problem. The numerical scheme uses polar slicing and a flux conservative version of 
the fluid equations of motion which are differenced using a two-step Lax-Wendroff 
scheme. The differencing scheme is second order accurate in both space and time; 
this has been verified by performing a sequence of evolutions at various levels of 
discretization. The code was used to confirm the results obtained by Evans and 
Coleman when k = 1/3; at the threshold of black hole formation self-similarity 
is observed over two orders of magnitude in scale, and the critical exponent is 
computed to be {3 ~ 0.352. 

For equations of state with < k <, 0.899, the critical exponents which have 
been computed from sequences of numerical evolutions agree well with those ob- 
tained using perturbation methods. Several values are tabulated in Table [j]. In 
near critical evolutions, continuous self-similarity is observed over approximately 
two orders of magnitude in scale before the matter either disperses or collapses into 
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Table 1: The scaling exponent for black hole mass for several values of k. The values predicted 
using perturbation theory are in the second column labelled /^predicted- The numerically observed 
values and estimates of their errors are also presented 



k 


/^predicted 


/^observed 


k 


/^predicted 


/^observed 


0.1 


0.1875 


0.196 ± 0.009 


0.6 


0.5556 


0.565 ± 0.011 


0.2 


0.2614 


0.265 ± 0.009 


0.7 


0.6392 


0.653 ± 0.005 


0.3 


0.3322 


0.337 ± 0.002 


0.8 


0.7294 


0.740 ± 0.003 


0.4 


0.4035 


0.408 ± 0.005 


0.9 




0.814 ± 0.001 


0.5 


0.4774 


0.486 ± 0.004 


0.95 




0.850 ± 0.005 



a black hole. 



1 

S 0.1 
0.01 



1 

J o.i 

0.01 



io- 3 10- 2 io- 1 10° io- 3 lo- 2 io- 1 10° 

i * I 

p - p 

Figure 2: The scaling law for black hole mass in marginally supercritical evolutions. Results of 
a best fit to a power law Mgjj oc \rj — r]*\P are presented for several values of k > 0.8. There is 
evidence for scaling over two orders of magnitude in black hole mass; the lower limit is set by the 
resolution of runs which had Ar = 0.001. 

The perturbative methods of Maison and Koike et aZ.i cannot be applied 
when k > 0.899 because globally analytic, self-similar solntkms representing perfect 
fluid collapse fail to exist for such stiff equations of state.tirEj Nevertheless, Type II 
phase transitions are observed in numerical simulations of collapse for these values 
of k. The code allows resolution of black holes down to ~ 10 _3 M, where M is the 
ADM mass on the initial slice. Figure || shows the scaling of black hole mass for 
several values of k > 0.8. 

Near criticality, the solution develops a highly evacuated region outside the 
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Figure 3: Profiles of f(r, t) when k = 0.95 showing the self-similarity over approximately two orders 
of magnitude in scale. The intermediate four profiles are overlayed when scaled in accordance with 

the self-similar ansatz. 

collapsing matter, however there is evidence for self-similarity over two orders of 
magnitude in scale in slightly sub-critical solutions. Several profiles of / [see Eq. 
(pi)] are shown in Fig ||. Some of the profiles are overlayed when scaled according 
to the self-similar ansatz; continuous self-similarity is apparent. Since stiff fluid 
solutions can be recast as scalar field solutions, one might expect that DSS should 
turn on as k — > 1. No evidence to this effect has been found in the numerical 
simulations. The reason appears to be connected with a lack of analyticity at the 
outer boundary of the collapsing material in perfect fluid collapse. 

2.5 Cosmic censorship and critical collapse 

The scaling law for black hole mass in Type II phase transitions suggests that black 
holes of arbitrarily small mass can be formed in gravitational collapse. At the crit- 
ical point, the central region of collapsing material is self-similar; an infinite train 
of scale invariant echoes accumulates within finite proper time in the echoing solu- 
tions. These two observations suggest that a singularity must form at the center p£ 
critical solutions. Is it naked? Numerical simulations by Hamade and Stewart lj 
provide direct evidence that it is. Using an evolution scheme based on double null 
coordinates, they find a regular Cauchy horizon with almost vanishing flux of scalar 
field across it. Gundlach has analytically continued the critical solution, constructed 
by his pseudo-spectral method, to the Cauchy horizon which is regular. Further- 
more, self-|Similar, perfect fluid solutions are known to have naked singularities at 
the origin. til Thus, suitably chosen, regular initial data can form a naked singularity 
in gravitational collapse; however, the data are not generic in the sense that they 
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belong to a sub-space of codimcnsion one in the space of all initial dataJa 

Weak cosmic censorship states that, in the evolution of generic initial conditions, 
singularities are always hidden behind event horizons. Since critical solutions violate 
the generic condition, they are not counter-examples to cosmic censorship. Two 
further assumptions arc implicit in the above statement: First, gravitational collapse 
is governed by the classical Einstein equations. Second, matter is described by 
fundamental fields on spacetime. Thus the formation of naked singularities in the 
collapse of dust is not usually considered a violation of cosmic censorship since 
dust only provides an effective description of matter, and it can become singular in 
fiat spacetime. Nevertheless, it is remarkable that naked singularities are so easy 
to form in such matter. Jhingan et alE3 have examined dust collapse using the 
Tolman-Bondi-Lemaitre metric, relating the formation of black holes and naked, 
shell-focusing singularities in such a collapse to the generic form of regular initial 
data. Such data are characterized by the density and velocity profiles of the matter 
on some initial time slice. Given a generic initial density profile, they have shown 
that there exists a corresponding velocity field which gives rise to a strong curvature, 
naked singularity in the evolution. This establishes that strong naked singularities 
arise for generic density profiles in the spherically symmetric collapse of dust. 

3 Quantum effects in critical spacetimes 

It is well known, since the work of Hawking, that a black hole decreases in size by 
emitting particles via quantum processes. The radiated particles have a thermal 
spectrum; for a static black hole the temperature is inversely proportional to the 
mass of the black hole. Thus, the smaller the black hole, the more it radiates. A 
black hole of mass M radiates all its mass in approximately 10~ 27 (A//lg) 3 seconds. 
Black holes formed in marginally super-critical collapse will therefore evaporate 
almost instantaneously; quantum effects are important in a full description of Type 
II phase transitions. 

In the absence of a complete theory of quantum gravity, semi-classical calcula- 
tions can provide some information about the underlying quantum evolution. Since 
renormalization breaks conformal invariance, scale invariant critical solutions surely 
are modified when curvatures approach Planck scales. Furthermore, the introduc- 
tion of a fundamental length scale, the Planck length, into the picture suggests that 
a mass-gap might occur at the threshold of black hole formation in a semi-classical 
treatment of critical point behavior. These speculations have been confirmed in 
several model problems. 

3.1 Two-dimensional dilaton models 

Peleg et aZ0 have studied one-loop quantum effects on the collapse of a massless 
scalar field in two-dimensional (2D) dilaton gravity. Reflecting boundary conditions 
are imposed at some finite value of the dilaton = <j) c in order to avoid the strong 
coupling regime of the theory which has a timelike singularity where exp(20) — > oo. 
The classical solutions exhibit critical point behavior: in supercritical evolutions 
the black hole mass -Mbh scales as Mbh oc \rj — r/»| /3 near to criticality. In this 2D 
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model, the critical exponent is (3 ~ 0.53. 

There is some freedom in constructing the effective action which describes the 
semi-classical theory since one is free to add local counter-terms to the Polyakov- 
Liouville term derived from the trace anomaly in two dimensions. Pclcg et pi add 
a term which makes the theory exactly soluble at the semi-classical level£J The 
quantum coupling constant k = N%/Y1 depends on the number of fields N. As 
expected, quantum effects are not relevant in the formation of sufficiently large 
black holes. The classical scaling relation holds over four orders of magnitude when 
k = 0.001, and over two orders of magnitude when k = 0.01. The threshold of black 
hole formation is characterized by a mass-gap in semi-classical evolutions, i.e. Mbh 
approaches a non-zero lower limit as 7/ tends to its critical value from above. The 
mass-gap depends both on k, and on the initial data. This is not surprising since 
the semi-classical action is non-local. Peleg has constructed analytic arguments 
supporting these results (see his contribution to this volume). 

3.2 Semi- classical models of spherical collapse 

The benefit of 2D models is that the quantum stress_energy tensor can be calcu- 
lated exactly with very little effort. Ayal and Piran c3 have considered the four- 
dimensional, semi-classical Einstein equations 

G = 8tt(T+(T)) (15) 

where G is the Einstein tensor, T is the stress energy of classical matter, and (T) is 
the renormalized stress-energy tensor of quantum matter. Since direct computation 
of (T) is extremely difficult, if not beyond current techniques, they use a model 
stress-energy tensor constructed from a two-dimensional one computed on the radial 
two sections d of the four dimensional spacetime. Tangential stresses are ignored, 
and (in a slight abuse of notation) 4nr 2 (T) = [1 + (a/r 2 ) 2 ]^ 1 (T( 2 )) where y/a 
is a length scale of order the Planck scale. The prefactor is required so that the 
stress-energy tensor is not singular at the origin in four dimensions. Even though 
this term leads to small violations of covariant conservation, it is a useful model of 
semi-classical effects including Hawking evaporation. 

Solving the semi-classical equations ([D]) numerically for scalar field collapse, 
Ayal and Piran show that quantum effects reduce the mass of a black hole compared 
to its classical value. When the length scale ^fa is increased the black hole no longer 
forms. As one might expect, the results are qualitatively similar to those of Peleg 
et al. 

3.3 Four- dimensional model calculation 

An alternative approach exploits the symmetry of spacetimes with continuous self- 
similarity to compute the general form, of the renormalized stress-energy tensor 
(RSET) for conformally coupled fields.EJ Spherically, symmetric self-similar space- 
times are conformally static, therefore one can apply the transformation rule for 

6y radial two section we mean the 2D spacetime obtained by ignoring the angular dimensions 
in the spherical line element. 
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the RSET, as derived by Page,@ to show that 

(T) = hr 4 [ functions of £ = r/t only] . (16) 

Here it is assumed that t is the coordinate in Eq. (|3|) normalized so that t = at 
the singularity in the critical solution. 



is 




Figure 4: The horizon location is determined by the roots of the function F(rh,£h) m Eq. jl^). 
We show here a schematic representation for several values of r\ which determine deviations from 
classically critical initial data, and A > 2. For sufficiently large r\ classical collapse takes hold and 
a black hole forms at R c . The function has a minimum, however, and another root R q exists. As 
rj is tuned to a critical value r\ q the two roots coincide. When r\ < r\ q no black hole forms. 



The critical solution in perfect fluid collapse is / = /*(£)■ The classical argu- 
ment which relates the Lyapunov exponent to the scaling relation for black hole 
mass can now be repeated including perturbations which originate from the quan- 
tum stress-energy tensor in Eq. (jl^). The semi-classical equations imply that the 
perturbations to the self-similar solution can be written as 

F(li,t) = m)-(r ] -r 1 *)fi(Ot- x +hf q (Ot- 2 , (17) 

where the last term arises from quantum effects. Clearly, when A < 2 the quantum 
perturbations will dominate as t — * 0. Thus, we restrict attention to equations of 
state with k <; 0.53 for which A < 2. (Sec Brady and Ottewill til for a discussion of 
k < 0.53.) Now, substitute t = r/£ into Eq. (|l7j). For fixed £ = £h the solutions of 

F(th,r h ) - /„(&) -(v- V*)h(Sh)(r h /S h )- x + ^fM)(rh/^h)- 2 = (18) 



determine the radii of apparent horizons in spacetime as a function of rj. Assuming 
that quantum effects compete with classical perturbations to reduce the mass of a 
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black hole that forms in gravitational collapse, F(^h,fh) is plotted schematically 
in Fig. ||. For 77 sufficiently greater than 77*, gravitational collapse and black hole 
formation is dominated by the classical perturbations; the apparent horizon is at 
R c . Note that the function has a minimum and a second root at R q < R c . As r\ 
decreases R c — > R q until the roots coincide at some r\ q . When r\ < r\ q no black hole 
forms; therefore, we infer that quantum effects induce a mass gap at the threshold of 
black hole formation in critical collapse of perfect fluids. It is also apparent that the 
classical scaling law for black-hole mass remains valid for black holes with apparent 
horizons significantly above the Planck length. 

4 Concluding remarks 

Critical point behavior in gravitational collapse is an exciting area of research in 
general relativity, one which continues to bring new insights into gravitational col- 
lapse and black hole formation. Significant progress has been made to understand 
the phenomenology observed in near critical evolutions. Nevertheless, important 
open questions remain. For example, is cosmic censorship upheld in all families of 
interpolating solutions for realistic matter fields? Or, when both Type I and Type 
II phase transitions occur, which properties of the initial data determine the critical 
point behavior? Numerical studies of gravitational collapse can provide a great deal 
of insight into the answers to these and other questions, however real progress will 
be made at the interface between numerical and mathematical approaches. 
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